All-optical mode division demultiplexing

ABSTRACT

Methods of decoding N superimposed coherent optical transmission modes transmitted along a multimode optical fiber are provided where the optical signal detector comprises coherent detection hardware and a digital signal processor. The optical signals are split into N optical detection channels of the coherent detection hardware, which is used to measure the phase and amplitude of the output optical signals as a vector field matrix [ê] 1×n . The digital signal processor determines the output principal states [PS out ] 1×n  of the output optical signals from the principal state eigenvectors of an output matrix corresponding to the N propagating optical signals at the receiving portion of the data transmission link and extracts the input optical signals [M] In  from the output optical signals using the output principal states [PS out ] 1×n , the vector field matrix [ê] 1×n  of the output optical signals, and a time delay τ of each output signal. Optical signal receivers for the execution of the methodology disclosed and contemplated herein are also provided.

This application claims the benefit of priority under 35 U.S.C. §119 of U.S. Provisional Application Ser. No. 62/030,182 filed on Jul. 29, 2014, the content of which is relied upon and incorporated herein by reference in its entirety.

BACKGROUND

Multimode optical fibers have been introduced to increase transmission capacity in optical communications networks. Wavelength division multiplexing is utilized in these multi-channel fibers to enable the transmission of multiple optical modes in a single core or multiple distinct cores and help meet the explosive growth in the demands on network capacity. However, mode coupling, which causes the signals carried by specific modes to couple among each other, is a challenging technical issue in this context. A variety of approaches have been introduced to deal with the information scrambling resulting from this mode coupling.

Multiple-Input Multiple-Output communication in multimode optical fiber communication channels (MIMO) relies on coherent transmission with extensive fast forward error correction and electronic compensation to sort out scrambled signals. The principal states approach, as set forth in US PG Pub. No. 2014/0161439 A1, seeks to launch light into the normal modes of the optical system to limit mode coupling. Unfortunately, the components needed for generating, combining, and splitting these states are relatively complex and can be challenging to update and control.

BRIEF SUMMARY

The present disclosure relates primarily to a system and methodology by which N output modes of an optical communications link can be deconvolved to extract the input communication signals. More particularly, the multiple wavelength signals from the N simultaneously received modal signals can be decoded by measuring the output principal states and the received modal fields rather than modal power information.

According to some embodiments of the present disclosure, methods of decoding N superimposed coherent optical transmission modes transmitted along a multimode optical fiber, are provided where the method comprises (i) encoding, combining, and transmitting N coherent input optical signals (M1, M2, M3, . . . ) along a multimode optical fiber at a launching portion of a data transmission link, wherein the launching portion of the data transmission link comprises a data encoder and a mode combiner; (ii) detecting each of the N propagating optical signals as output optical signals at a receiving portion of the data transmission link, wherein the receiving portion of the data transmission link comprises coherent detection hardware and a digital signal processor; (iii) splitting the output optical signals into N optical detection channels of the coherent detection hardware; (iv) utilizing the coherent detection hardware to measure the phase and amplitude of the output optical signals as a vector field matrix [ê]_(1×n) associated with the output modes, M; (v) utilizing the digital signal processor to determine output principal states [PS_(Out)]_(1×n) of the output optical signals from the principal state eigenvectors of an output matrix corresponding to the N propagating optical signals at the receiving portion of the data transmission link; and (vi) utilizing the digital signal processor to extract the input optical signals [M]_(In) from the output optical signals using the output principal states [PS_(Out)]_(1×n), the vector field matrix [ê]_(1×n) of the output optical signals, and a time delay τ of each output signal.

According to other embodiments of the present disclosure, optical signal receivers are provided for the execution of the methodology disclosed and contemplated herein.

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS

The following detailed description of specific embodiments of the present disclosure can be best understood when read in conjunction with FIG. 1, which is a schematic illustration of a data transmission link that may be implemented in accordance with the concepts of the present disclosure.

DETAILED DESCRIPTION

FIG. 1 is a schematic illustration of a data transmission link 100 that may be implemented in accordance with the concepts of the present disclosure. It is contemplated that the data transmission link 100 may be configured as all or part of data transmission or telecommunication network, with varying degrees of complexity. Generally, in FIG. 1, respective coherent transmission modes M1, M2, M3, etc., are encoded with data and input at a launching portion 10 of the link 100 using a data encoder 20, which allows attributes of the respective coherent transmission modes M1, M2, M3, etc., to be controlled separately, and a mode combiner 30, which multiplexes the encoded transmission modes M1, M2, M3, etc., for propagation along a multimode optical fiber 40 via suitable fiber links 45. Although the data encoder 20 is illustrated schematically in FIG. 1 with reference to separate control of the phase and amplitude of N superimposed coherent transmission modes M1, M2, M3, etc., it is contemplated that the data encoder 20 may be additionally or alternatively configured to encode the polarization and/or other attributes of the N superimposed coherent transmission modes M1, M2, M3, etc. It is contemplated that the launching portion 10 of the data transmission link 100 of the present disclosure may be utilized to implement a variety of encoding schemes including, for example, a coherent optical MIMO (multiple-input, multiple-output) system.

The propagating signals are subsequently received and decoded at a receiving portion 50 of the data transmission link 100, which may also be referred to as an optical signal receiver and is illustrated schematically to comprise coherent detection hardware 60 (including a local oscillator LO) and a digital signal processor (DSP) 70. Analog-to-digital converters A/D are also provided to enable digital signal processing. The particular configuration of the coherent detection hardware 60 is beyond the scope of the present disclosure and may be gleaned from any of a variety of conventional or yet-to-be developed teachings on the subjects. See, for example, Ly-Gagnon et al., “Coherent Detection of Optical Quadrature Phase-Shift Keying Signals with Carrier Phase Estimation,” Journal of Lightwave Technology, Vol. 24, No. 1, pp. 12-21, January 2006. The DSP 70 is programmed to cooperate with the coherent detection hardware 60 and to execute the signal processing algorithms of the present disclosure. Beyond this, the particular configuration of the DSP 70 is also beyond the scope of the present disclosure. Zhang, “Digital Signal Processing for Optical Coherent Communication Systems,” DTU Fotonik, 2012, presents useful teachings on the subject that complement the present disclosure. Particular reference is made to Chapter 3 of Zhang: “Coherent Detection with DSP Algorithms.”

US PG Pub. No. 2014/0161439 A1, teaches suitable hardware and methodology for generating principal state eigenvectors of an output matrix corresponding to N propagating optical signals at a receiving portion of a data transmission link. The aforementioned patent publication also teaches a matrix method that can be used to solve for the eigenvalues and eigenvectors of the principal states of the N propagating optical signals. These principal sate eigenvalues and eigenvectors can be used to facilitate the signal decoding described herein.

More specifically, the present disclosure presents methodology and associated architecture for decoding N superimposed coherent optical transmission modes transmitted along a multimode optical fiber, in which mode coupling occurs. Referring to FIG. 1, given N coherent input optical signals (M1, M2, M3, . . . ) encoded, combined, and transmitted along a multimode optical fiber 40 at the launching portion 10 of the data transmission link 100, where N>1, each of the N propagating optical signals are detected as output optical signals at the receiving portion 50 of the data transmission link 100 by splitting the output optical signals into N optical detection channels of the coherent detection hardware 60 of the receiving portion 50.

The coherent detection hardware 60 is used to measure the phase and amplitude of the output optical signals as a vector field at each modal position designated with a matrix [ê]_(1×n) and the digital signal processor 70 determines the output principal states [PS_(Out)]_(1×n) of the output optical signals from the principal state eigenvectors of an output matrix corresponding to the N propagating optical signals at the receiving portion 50 of the data transmission link 100. The digital signal processor 70 subsequently extracts the input optical signals [M]_(In) from the output optical signals using the output principal states [PS_(Out)]_(1×n), the vector field matrix [ê]_(1×n) of the output optical signals, and a time delay τ of each output signal. This extraction process is described from the three different perspectives herein: first, with a general reference to calculating the complex transpose of the output principal states of polarization (PSPs); second, with general reference to deconvolving the input principal states; and third, with a more specific reference to using [γ]_(n×n).

More specifically, the digital signal processor can be used to determine the superimposed input states [S_(In)]_(1×n) of the input optical signals [M]_(In) from the output principal states [PS_(Out)]_(1×n) of the output optical signals [M]_(Out) by calculating the complex transpose of the output principal states [PS_(Out)]_(1×n). The generalized transfer matrix [H] of the system can then be calculated from the superimposed input states [S_(In)]_(1×n) and the output principal states [PS_(Out)]_(1×n) according to:

[PS _(Out)]_(1×n) =[H][S _(In)]_(1×n).

M1, M2, M3, . . . can then be extracted based at least in part on the relation

[M] _(In) =[H′] ⁻¹ [M] _(Out)

where the input optical signals [M]_(In) are represented as follows:

[M] _(In)=(M1, M2, M3, . . . ).

From a different perspective, the methodology and architecture of the present disclosure may be described with reference to the use of the digital signal processor 70 to initially determine the superimposed input states [S_(In)]_(1×n) of the input optical signals [M]_(In) from the output principal states [PS_(Out)]_(1×n) of the output optical signals [M]_(Out). Subsequently, the input optical signals [M]_(In) are determined by deconvolving [S_(In)]_(1×n) at least partially according to

[M] _(In) =[D][S _(In)]_(1×n); and

M1, M2, M3, . . . are extracted from [M]_(In), where [M]_(In)=(M1, M2, M3, . . . ), and [D] is a matrix establishing a relationship between the superimposed input states and the combined input optical signals

From another perspective, the methodology and architecture of the present disclosure may be described with reference to the initial determination of [a]_(n×n) from the measured vector field matrix [ê]_(1×n) and the output principal states [PS_(Out)]_(1×n) at least partially according to

[PS _(Out)]_(1×n) =[a] _(n×n) [ê] _(1×n).

Where [a]_(n×n) relates the output modes to the output principal states as described in US2014/016,439 A1, hereby incorporated by reference in its entirety. The superimposed input states [S_(In)]_(1×n) can be subsequently related to the vector fields f_(In) of the input optical signals (M1, M2, M3, . . . ) at least partially according to

[f _(In) S _(In)]_(1×n)=[γ]_(n×n) [ê] _(1×n)

where the elements of [γ]_(n×n) are constructed from [a]_(n×n), time dependent vector fields f_(j) of the output optical signals, and propagation constants β_(i) associated with the time dependent vector fields f_(j). Accordingly, M1, M2, M3, . . . can be determined at least partially from the vector fields f_(In), where [M]_(In)=(M1, M2, M3, . . . ). From this perspective, the digital signal processor 70 can be programmed to determine the elements of [γ]_(n×n) at least partially according to

$\lbrack\gamma\rbrack_{nxn} = \begin{bmatrix} {a_{11}e^{{- j}\; \beta_{1}z}f_{1}} & {a_{12}e^{{- j}\; \beta_{1}z}f_{2}} & \ldots \\ {a_{21}e^{{- j}\; \beta_{2}z}f_{1}} & {a_{22}e^{{- j}\; \beta_{2}z}f_{2}} & \vdots \\ \vdots & \ldots & \ddots \end{bmatrix}$

where a₁₁, a₁₂, . . . are elements of [a]_(n×n), j is the imaginary number, β₁, β₂, . . . , are the propagation constants associated with the principal states time delays, i.e., the eigenvalues of the principal states eigenvalue equation, z is the length of the optical fiber 40, and f₁, f₂, . . . are the time dependent vector fields of the output optical signals. In one embodiment, it is contemplated that the digital signal processor 70 can be programmed to determine the propagation constants β_(i) at least partially according to

β_(i)=τ_(i) kc

where k=2π/λ and c is the speed of light.

To determine the input modulated fields, the equation

[f _(In) S _(In)]_(1×n)=[γ]_(n×n) [ê] _(1×n)

can be expanded to give

S1_(In) =a ₁₁ e ^(−jβ) ¹ ^(z) f ₁ {circumflex over (m)} ₁ +a ₁₂ e ^(−jβ) ¹ ^(z) f ₂ {circumflex over (m)} ₂ +a ₁₃ e ^(−jβ) ¹ ^(z) f ₃ {circumflex over (m)} ₃+ . . .

S2_(In) =a ₂₁ e ^(−jβ) ² ^(z) f ₁ {circumflex over (m)} ₁ +a ₂₂ e ^(−jβ) ² ^(z) f ₂ {circumflex over (m)} ₂ +a ₂₃ e ^(−jβ) ² ^(z) f ₃ {circumflex over (m)} ₃+ . . .

S3_(In) =a ₃₁ e ^(−jβ) ³ ^(z) f ₁ {circumflex over (m)} ₁ +a ₃₂ e ^(−jβ) ³ ^(z) f ₂ {circumflex over (m)} ₂ +a ₃₃ e ^(−jβ) ³ ^(z) f ₃ {circumflex over (m)} ₃+ . . .

where {circumflex over (m)} represents the specific polarization mode to be detected. The resulting matrix can be diagonalized such that the diagonal components of the diagonalized matrix are the input modulated fields.

As a result of this method, by using this algorithm and the measured output principal states as in US 2014/016,439 A1, one can deconvolve N output signals using time independent mode multiplexing and demultiplexing components, thereby eliminating the need to provide time dependent multiplexing and demultiplexing as described in US 2014/016,439 A1. These components, as described therein, are required to multiplex and demultiplex the superimposed principal states which will vary in time.

It is noted that recitations herein of a component of the present disclosure being “programmed” or “configured” in a particular way, to embody a particular property, or function in a particular manner, are structural recitations, as opposed to recitations of intended use. More specifically, the references herein to the manner in which a component is “programmed” or “configured” denotes an existing physical condition of the component and, as such, is to be taken as a definite recitation of the structural characteristics of the component.

Having described the subject matter of the present disclosure in detail and by reference to specific embodiments thereof, it is noted that the various details disclosed herein should not be taken to imply that these details relate to elements that are essential components of the various embodiments described herein, even in cases where a particular element is illustrated in each of the drawings that accompany the present description. Rather, the claims appended hereto should be taken as the sole representation of the breadth of the present disclosure and the corresponding scope of the various inventions described herein. Further, it will be apparent that modifications and variations are possible without departing from the scope of the invention defined in the appended claims. More specifically, although some aspects of the present disclosure are identified herein as preferred or particularly advantageous, it is contemplated that the present disclosure is not necessarily limited to these aspects. 

1. A method of decoding N combined coherent optical transmission modes transmitted along a multimode optical fiber, where N>1, and the method comprises: encoding, superimposing, and transmitting N coherent input optical signals (M1, M2, M3, . . . ) along a multimode optical fiber at a launching portion of a data transmission link, wherein the launching portion of the data transmission link comprises a data encoder and a mode combiner; detecting each of the N propagating optical signals as output optical signals at a receiving portion of the data transmission link, wherein the receiving portion of the data transmission link comprises coherent detection hardware and a digital signal processor; splitting the output optical signals into N optical detection channels of the coherent detection hardware; utilizing the coherent detection hardware to measure the phase and amplitude of the output optical signals as a vector field matrix [ê]_(1×n); utilizing the digital signal processor to determine output principal states [PS_(Out)]_(1×n) of the output optical signals from the principal state eigenvectors of an output matrix corresponding to the N propagating optical signals at the receiving portion of the data transmission link; and utilizing the digital signal processor to extract the input optical signals [M]₁, from the output optical signals using the output principal states [PS_(Out)]_(1×n), the vector field matrix [ê]_(1×n) of the output optical signals, and a time delay τ of each output signal.
 2. A method as claimed in claim 1 wherein the method further comprises: utilizing the digital signal processor to determine superimposed input states [S_(In)]_(1×n) of the input optical signals [M]_(In) from the output principal states [PS_(Out)]_(1×n) of the output optical signals [M]_(Out) by calculating the complex transpose of the output principal states [PS_(Out)]_(1×n); utilizing the digital signal processor to calculate the transfer matrix [H] of the system from the superimposed input states [S_(In)]_(1×n) and the output principal states [PS_(Out)]_(1×n), where [PS _(Out)]_(1×n) =[H][S _(In)]_(1×n), and to extract M1, M2, M3, . . . based at least in part on the relation [M] _(In) =[H′] ⁻¹ [M] _(Out) where [M] _(In)=(M1, M2, M3, . . . ).
 3. A method as claimed in claim 1 wherein the method further comprises utilizing the digital signal processor to: determine superimposed input states [S_(In)]_(1×n) of the input optical signals [M]_(In) from the output principal states [PS_(Out)]_(1×n) of the output optical signals [M]_(Out); determine the input optical signals [M]_(In) by deconvolving [S_(In)]_(1×n) at least partially according to [M] _(In) =D[S _(In)]_(1×n); and extract M1, M2, M3, . . . from [M]_(In), where [M]_(In)=(M1, M2, M3, . . . ).
 4. A method as claimed in claim 1 wherein the digital signal processor: determines [a]_(n×n) from the measured vector field matrix [ê]_(1×n) and the determined output principal states [PS_(Out)]_(1×n) at least partially according to [PS _(Out)]_(1×n) =[a] _(n×n) [ê] _(1×n); relates the superimposed input states [S_(In)]_(1×n) to vector fields f_(In) of the input optical signals (M1, M2, M3, . . . ) at least partially according to [f _(In) S _(In)]_(1×n)=[γ]_(n×n) [ê] _(1×n) where the elements of [γ]_(n×n) are constructed from [a]_(n×n), time dependent vector fields f_(j) of the output optical signals, and propagation constants β_(i) associated with the time dependent vector fields f_(j); and determines M1, M2, M3, . . . at least partially from the vector fields f_(In), where [M]_(In)=(M1, M2, M3, . . . ).
 5. A method as claimed in claim 4 wherein the digital signal processor determines the elements of [γ]_(n×n) at least partially according to $\lbrack\gamma\rbrack_{nxn} = \begin{bmatrix} {a_{11}e^{{- j}\; \beta_{1}}f_{1}} & {a_{12}e^{{- j}\; \beta_{1}}f_{2}} & \ldots \\ {a_{21}e^{{- j}\; \beta_{2}}f_{1}} & {a_{22}e^{{- j}\; \beta_{2}}f_{2}} & \vdots \\ \vdots & \ldots & \ddots \end{bmatrix}$ where a₁₁, a₁₂, . . . are elements of [a]_(n×n), β₁, β₂, . . . , are the propagation constants associated with the principal states time delays, and f₁, f₂, . . . are the time dependent vector fields of the output optical signals.
 6. A method as claimed in claim 5 wherein the digital signal processor determines the input optical signals (M1, M2, M3, . . . ) by: expanding [f _(In) S _(In)]_(1×n)=[γ]_(n×n) [ê] _(1×n) to give S1_(In) =a ₁₁ e ^(−jβ) ¹ ^(z) f ₁ {circumflex over (m)} ₁ +a ₁₂ e ^(−jβ) ¹ ^(z) f ₂ {circumflex over (m)} ₂ +a ₁₃ e ^(−jβ) ¹ ^(z) f ₃ {circumflex over (m)} ₃+ . . . S2_(In) =a ₂₁ e ^(−jβ) ² ^(z) f ₁ {circumflex over (m)} ₁ +a ₂₂ e ^(−jβ) ² ^(z) f ₂ {circumflex over (m)} ₂ +a ₂₃ e ^(−jβ) ² ^(z) f ₃ {circumflex over (m)} ₃+ . . . S3_(In) =a ₃₁ e ^(−jβ) ³ ^(z) f ₁ {circumflex over (m)} ₁ +a ₃₂ e ^(−jβ) ³ ^(z) f ₂ {circumflex over (m)} ₂ +a ₃₃ e ^(−jβ) ³ ^(z) f ₃ {circumflex over (m)} ₃+ . . . where {circumflex over (m)} represents the specific polarization mode to be detected; and diagonalizing the resulting matrix such that the diagonal components of the diagonalized matrix represent input modulated fields of the input optical signals (M1, M2, M3, . . . ).
 7. A method as claimed in claim 4 wherein the digital signal processor determines the propagation constants β_(i) at least partially according to β_(i)=τkc where k=2π/λ and c is the speed of light.
 8. An optical signal receiver comprising coherent detection hardware and a digital signal processor programmed to cooperate with the coherent detection hardware to decode N superimposed coherent optical transmission modes transmitted along a multimode optical fiber, where N>1, wherein the decoding comprises: detecting each of the N propagating optical signals as output optical signals; splitting the output optical signals into N optical detection channels of the coherent detection hardware; utilizing the coherent detection hardware to measure the phase and amplitude of the output optical signals as a vector field matrix [ê]_(1×n); utilizing the digital signal processor to determine output principal states [PS_(Out)]_(1×n) of the output optical signals from the principal state eigenvectors of an output matrix corresponding to the N propagating optical signals at the receiving portion of the data transmission link; and utilizing the digital signal processor to extract the input optical signals [M]_(In) from the output optical signals using the output principal states [PS_(Out)]_(1×n), the vector field matrix [ê]_(1×n) of the output optical signals, and a time delay τ of each output signal.
 9. An optical signal receiver as claimed in claim 8 wherein the digital signal processor is programmed such that the decoding further comprises: utilizing the digital signal processor to determine superimposed input states [S_(In),]_(1×n) of the input optical signals [M]_(In) from the output principal states [PS_(Out)]_(1×n) of the output optical signals [M]_(Out) out by calculating the complex transpose of the output principal states [PS_(Out)]_(1×n); utilizing the digital signal processor to calculate the transfer matrix [H] of the system from the superimposed input states [S_(In)]_(1×n) and the output principal states [PS_(Out)]_(1×n), where [PS _(Out)]_(1×n) =[H][S _(In)]_(1×n), and to extract M1, M2, M3, . . . based at least in part on the relation [M] _(In) =[H′] ⁻¹ [M] _(Out) where [M] _(In)=(M1, M2, M3, . . . ).
 10. An optical signal receiver as claimed in claim 8 wherein the digital signal processor is programmed such that the decoding further comprises: determining superimposed input states [S_(In)]_(1×n) of the input optical signals [M]_(In) from the output principal states [PS_(Out)]_(1×n) of the output optical signals [M]_(Out); determining the input optical signals [M]_(In) by deconvolving [S_(In)]_(1'n) at least partially according to [M] _(In) =D[S _(In)]_(1×n); and extracting M1, M2, M3, . . . from [M]_(In), where [M]_(In)=(M1, M2, M3, . . . ).
 11. An optical signal receiver as claimed in claim 8 wherein the digital signal processor is programmed such that the decoding further comprises: determining [a]_(n×n) from the measured vector field matrix [ê]_(1×n) and the determined output principal states [PS_(Out)]_(1×n) at least partially according to [PS _(Out)]_(1×n) =[a] _(n×n) [ê] _(1×n); relating the superimposed input states [S_(In)]_(1×n) to vector fields f_(In) of the input optical signals (M1, M2, M3, . . . ) at least partially according to [f _(In) S _(In)]_(1×n)=[γ]_(n×n) [ê] _(1×n) where the elements of [γ]_(n×n) are constructed from [a]_(n×n,) time dependent vector fields f_(j) of the output optical signals, and propagation constants β_(i) associated with the time dependent vector fields f_(j); and determining M1, M2, M3, . . . at least partially from the vector fields f_(In), where [M]_(In)=(M1, M2, M3, . . . ).
 12. An optical signal receiver as claimed in claim 11 wherein the digital signal processor determines the elements of [γ]_(n×n) at least partially according to $\lbrack\gamma\rbrack_{nxn} = \begin{bmatrix} {a_{11}e^{{- j}\; \beta_{1}}f_{1}} & {a_{12}e^{{- j}\; \beta_{1}}f_{2}} & \ldots \\ {a_{21}e^{{- j}\; \beta_{2}}f_{1}} & {a_{22}e^{{- j}\; \beta_{2}}f_{2}} & \vdots \\ \vdots & \ldots & \ddots \end{bmatrix}$ where a₁₁, a₁₂, . . . are elements of [a]_(n×n), b₁, b₂, . . . , are the propagation constants associated with the time dependent vector fields of the output optical signals, and f₁, f₂, . . . are the time dependent vector fields of the output optical signals.
 13. An optical signal receiver gas claimed in claim 11 wherein the digital signal processor determines the propagation constants β_(i) at least partially according to β_(i)=τ_(i)kc where k=2π/λ and c is the speed of light.
 14. An optical communications system using the method as in claim 1, and the measured output principal states as described in US 2014/016,439 A 1 to deconvolve N output signals using mode multiplexing and demultiplexing components 